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ELASTIC INSTABILITY OF FLAT PLATES SUBJECTED 
TO PARTIAL EDGE LOADs! 


Lev Zetlin,? A.M. ASCE 


SYNOPSIS 


Buckling loads are determined of flat rectangular plates subject to loads 
distributed over short portions of the edges. The equations governing the 
buckling loads, which are partial differential equations with variable coef- 
ficients, are discussed and an approach to the solution is made by means of 
the Energy Method. Results of calculations have been reduced to curves 
suitable for design purposes. These curves are given for a wide range of 
plate dimensions and loading widths. The theoretical results are compared 
with a number of tests. 


Notation 


- the horizontal length (span) of plate. 
- the vertical height (depth) of plate. 
- h/L; non-dimensional ratio. 

~ thickness of plate. 


- width through which the external load is applied to the top edge 
of plate. 


- stress in the x-direction. 

- stress in the y-direction. 

- shearing stress 

- yield stress in the steel of the test specimen. 
An, Bn, Cp, and Dy, = - coefficients in the stress function. 


Amk OF Apr - (or A with any numerical subscript), amplitudes of 
the double Fourier series in the eigenfunction. 


D - flexural rigidity. 


w(x,y) - displacement of a point anywhere on the neutral plane of the plate 
in the direction perpendicular to the plane of the plate. 


. Based on a thesis by L. Zetlin submitted to Cornell University in June 
1953 in partial fulfillment of the requirements for the degree of Doctor of 
Philosophy. 

2. Asst. Prof., Structural Engineering, Cornell University, Ithaca, N. Y. 
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INTRODUCTION 


This paper presents an analysis of the elastic stability of a plate supported 
at both ends and loaded over a portion of one of the edges. The plate is shown 
in Fig. 1. Its edges are simply supported. The total edge load, P, is applied 
along the top edge and is uniformly distributed over a width 2B with intensity 
p = P/2B. The reactions are supplied by shear stresses with parabolic dis- 
tribution along the two vertical edges. Both the reactions and the external 
load P are symmetrical about the center line of the plate. The solution is for 
the critical value of p at which the plate buckles and is in terms of the three 
parameters, 8 = h/L, 2B and D. Thus, the solution applies to a wide range of 
variations in the dimensions of the plate and width of the strip loading. This 
small deflection analysis is based on the usual assumptions of elasticity, 
homogeneity and isotropy. 

Since the mathematical form of the solution is too complex for direct use 
by a designer, a number of graphs have been drawn from which the buckling 
loads can be easily read. 

To the author’s best knowledge, the only investigation of a somewhat simi- 
lar problem is that by K. Girkmann.(23,4) His results however, apply only to 
plates with h/L 2 0.9, because of certain mathematical limitations in his initial 
assumptions. This will be pointed out more explicitly later in connection with 
the mathematical development of the problem. Moreover, Professor Girk- 
mann’s papers give the solution only in the form of a determinant, which has 
to be evaluated for any particular case. 

Only the most essential mathematical steps and results, which are neces- 
sary for the understanding of the theoretical procedure in the solution of the 
problem, are given in this paper. For the complete mathematical develop- 
ment of the solution the reader is referred to footnote 1. 


I. Theoretical Solution 


The elastic instability of any flat plate is governed by the eigenvalues of 
the following differential equation (see for instance Ref. 1, page 324): 


The coefficients of the derivatives on the right-hand side of Eqn. (1) are 
stresses throughout the plate due to the external loads; these stresses are in 
the plane of the plate and are assumed to be constant through the thickness of 
the plate. Obviously, for the case shown in Fig. 1 the stresses vary from 


795-2 


4 4 4 

Dw aw aw 

* Dy? 

/ ra ~ 


point to point in the plane of the plate and the differential equation (1) is one 
with variabie coefficients. 

In order to be able to cope with Eqn. (1), the stress distribution throughout 
the plate has to be determined first. Since a solution in terms of plate dimen- 
sions and width of the strip loading as parameters is desired, the stresses 
have to be given as functions of the coordinates (xy). 


a) Stress Distribution in the Plate. 

The stress distribution in the plate of Fig. 1 due to the external loads as 
shown, can be determined by the following stress function taken from Ref. (2) 
and (5): 


where Ay, Bn, Cp and Dy, are constants which depend on the boundary forces 
along the horizontal edges of the plate. The shear stress distribution along 
the vertical edges resulting from this stress function is close to parabolic. 
Along the top edge, the external load p = P/2B is expanded in a Fourier 
series. 


ply 


The boundary conditions on the stresses are: 
aty=0;7,,=0.. +... . Eqn, (4a) and by - - renee Eqn.(4b) 
aty=h;Txy=0.........Eqn. (4c) and by . Eqn.(4d) 


Differentiation of Eqn. (2) as is indicated below, together with the proper 
substitution of the boundary stresses of Eqns (4), yield the following expres- 
sions for the stresses in the plate. 


oo ~ 
Be 2 
n= 135... 
_ (2) 
n L L 
N=1,35 ... 
n=135.... Egn.(3) 
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h=435.... 
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Eqn. (56) 
n=435... 
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The constants A,, Bn, Cy and Dp in Eqn. (5) (which are the same constants 
as in Eqn. (2)) aregiven by Eqn. (6) below: 


~ n ‘ 
A, => Ki, [14 e 


C, = Pak, [Ot 


Egu. (6) 


/ 


7 Eau (7) 


b) Determination of the Buckling Load. 

Having expressions for the stresses in Eqns. (5), one way to determine the 
critical load p is to substitute Eqns. (5) into Eqn. (1); then, one of the eigen- 
values of Eqn. (1) would yield the critical load, probably in the form of an 
infinite series. The mathematical background on which such a procedure 
should be based, as well as a discussion of a mathematical proof that a solu- 
tion to Eqn. (1) exists, is treated at some length in footnote 1. Such a proce- 
dure, however, if at all practically possible, would be an ambitious undertak- 
ing. 

Another, more practical way, to determine the critical load p, is the use 
of energy relationships of the plate at buckling, and subsequent application of 
the Rayleigh-Ritz method. This procedure, too, is discussed extensively in 
footnote 1, and its explanation is beyond the scope of this paper. Only the 
essential steps and equations of this procedure will be pointed out below. 

If at the instant that the plate passes from stable to unstable equilibrium, 
A U denotes the increment of the total potential energy of the plate and the 
external loads, AV denotes the increment of the strain energy of the plate, 
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and AT denotes the increment of the potential energy of the external loads, 
then: 


4U=AV+AT +++... Eqn.(8) 


Expressing AV and AT in terms of plate dimensions, its material proper- 
ties, and the stresses in the plate, the following relationship results: 


(9) 


An expression similar to Eqn. (9) can be found on page 325 of Ref. (1) and 
its derivation on pages (311-314 of the same reference. 

Stresses bx, by and ¢,. in Eqn. (9) are, of course, given by Eqns. (5). In 
order to be able to 3 proceed with the solution of the problem, an expression 
for w as a function of x and y, i.e. the deflected shape of the plate at buckling, 
should be known or properly assumed. The closer the assumption is to the 
actual buckled shape of the plate, the closer will be the computed critical load 
p to the actual buckling load. 

In this respect, the author had the advantage of tests on a related problem 
before the theoretical solution was attempted. From these tests, observations 
of the shape of the plate at failure had been made. The tests will be discussed 
in more detail later in the paper. 

The available test observations indicated qualitatively that the shape of 
such plates at failure, for a large range of variations in the parameters 2B, 
and h/L, could be represented with good accuracy by the following trigonome- 
tric series: 


where A_, are constants to be determined by the Rayleigh-Ritz method. 
That Rs assumption in Eqn. (10) is satisfactory, will be demonstrated 


* It will be noticed that Eqn. (10) satisfies the boundary conditions of the 
plate in Fig. 1. 
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later in comparing computed critical loads with the failure loads observed in 
the tests. 

At this stage the reader will be referred to the statement made in the 
Introduction in connection with the limitation of Professor Girkmann’s work. 
Professor Girkmann(2) assumed the following expression for w: 


sm p+ Bsn 
SEY Fe TX 
+ (sin | os ( Z ) 


Eqn. (11) has more limited possibilities to represent the shape of a buckled 
plate than Eqn. (10). Hence, the critical loads yielded by Eqn. (11) are likely 
to be less accurate than those yielded by Eqn. (10). This will be particularly 
so with smaller values of h/L (i.e. longer plates) and smaller values of 2B 
(width of the strip loading). 

The next step is to evaluate the increment of the total potential energy AU 
in Eqn. (9). For this purpose, Eqn. (10) is differentiated as indicated by Eqn. 
(9); the differentiated expressions and the expressions for the stresses as 
given by Eqns. (5) are substituted into Eqn. (9), and, finally, integration of 
Eqn. (9) is performed. The process is quite lengthy and its details will not 
be given here. It will be found in footnote 1. The described operation gives 
an expression for AU in terms of the parameters 2B, h/L and D, and the six 
unknown constants 

Since the total potential energy at elastic instability has a stationary value, 
(A U)/OAm, = 0. Hence, in the final step of the solution, 4 U is differenti- 
ated with respect to each of the six constants A,,, and each derivative is 
equated to zero. Thus, six homogeneous algebraic equations in the six un- 
knowns A,,,, (i.e. Ayy, A13, A415, Ao1, A23, Ags) are obtained. These equa- 
tions have the following form: 


AU) _ Is 
7 
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In Eqn. (12a), 
p is the external load, 
dh D/L? 


ng as follows: 
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are expressions in the form of infinite series which 
enh bath on B nd ba. The expressions for I are quite lengthy; to give an 
idea of their form, J’; is reproduced below: 


n= 1,35... 
(n*-2 “at 


p2 435... 


nr (n= (772) ii Pesnis+4) 
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(24) 
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It will be noted that the unknowns in Eqn. (12a) are the constants A, and 


the external load p. 
The information that equations (12a) can supply is: 


1) Values of p for which the set of the six equations in (12a) has non zero 
solutions for the A’s. There will be, in general, six such distinct values 
of p, the smallest of which is the first buckling load. 

2) Ratio between the six constants A,,,. If this ratio is known for any 
particular plate (i.e. certain values of 2B, h/L and D), the shape of the 
buckled plate can be represented by substituting the ratios of A,, into 
Eqn. (10). Of course, as is the case in all instability problems, the 
actual deflections w will remain indeterminate as to magnitude. 


Of the above two items, the first one is of practical interest. The second 
item may be of interest in research work as a check on the initial assumption 
of the buckled surface (Eqn. (10)). 

In solving equations (12a), the smallest value of p has to be determined 
first; the corresponding ratios between A, can then be determined. For the 
set of the six equations (12a) to have solutions, the determinant of these equa- 
tions has to be equal to zero, or: 


G 

2 

Sp, As + , be 


The largest eigenvalue of the determinant in Eqn. (13) will yield the small- 
est buckling load pe. Denoting this eigenvalue by 72k, 


4 
m*D 
mk = Pep = jews 
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the following formula for the critical load is obtained: 


2 
Por =k 7 P 
2 
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Formula (15) has the same form as the familiar formula for the critical 
load of a plate uniformly loaded on the top and the bottom edges (see Ref. (1), 
page 329, Eqn. 213). The constant k in Eqn. (15), is, of course, different from 
the equivalent constant in that relatively elementary case and depends both on 
the dimensions of the plate (h/L) and the ratio of the width 2B of the strip 
loading to L (but not on the thickness of the plate). 


Il. Numerical Evaluation and Results. 


As is seen from Egn. (15), in order to obtain the critical load for a plate 
as shown in Fig. 1, the value of the constant k corresponding to the dimen- 
sions of the plate and the width of the strip loading must be known. The work 
connected with the evaluation of this constant in each case by conventional 
means is prohibitive. However, it would be possible to carry out appropriate 
computations on electronic computers and subsequently to plot the computed 
values of k in the form of graphs which may be readily useable for design 
purposes. 

Only nine values of k for different combinations of 8 and 2B, have been 
computed on ordinary desk computers. Before presenting the results, how- 
ever, a brief discussion on the convergence of the series in Eqn. (12e) as 
well as on the method of computing the largest eigenvalue of the matrix in 
Eqn. (13) will be made. 


a) Convergence of the Series. 

The elements of the matrix in Eqn. (13) are each given by an infinite series 
similar to that in Eqn. (12e). The series for each of the 21 elements of the 
matrix are different in form, but the summation in all of them is carried on 
n, as in Eqn. (12e). The series and their convergence are more fully dis- 
cussed in footnote 1. The following are the results of that discussion: 


1) All 21 expressions for J are alternating infinite series. 

2) Convergence is faster for larger values of both 8 and 2B, the influence 
of 8 on the convergence being more pronounced than that of 2B. 

3) The convergence for numerical evaluation is generally fair. For ex- 
ample, for the cases computed in this paper, the number of terms 
necessary to obtain an accuracy of 0.1% in the summation of the series 
were as follows: 

a) B = 2, one to four terms, 

b) B = 1, one to six terms, 

c) B = 1/4, two to eight terms. 

Variation in the number of terms corresponding to each value of @ as 
given above, depended on the value of 2B and on the type of the series, 
i.e. on the subscript of T. 


b) Evaluation of the Eigenvalues. 

Numerical values of the largest eigenvalue of such matrices can be com- 
puted with relative ease by the iteration method.(6) For the cases computed 
in this paper, four to six iterations were necessary to achieve an accuracy of 
1% in the eigenvalues. 


795-11 


| 
| 


c) Results of Numerical Computations. 
The nine cases for which k has been computed, are shown in Table I: 


Table I 


Computed 


Case value of 

No. 2B/L k 
1 2 1/65 378.00 
2 2 1/20 116.30 
3 2 1 10.67 
4 1 1/65 216.00 
5 1 1/20 67.30 
6 1 1 6.20 
7 1/4 1/65 559.00 
8 1/4 1/20 169.20 
9 1/4 1 26.60* 


As was already mentioned, in order to obtain one value of k it is required 
to compute 21 infinite series—one for each element of the symmetric matrix— 
and then to evaluate the largest eigenvalue of the resulting numerical matrix. 
The summation of the infinite series for each element of the matrices for the 
nine cases shown in Table I was carried out until the remainder of the series 
in each case was less than 0.1% of the sum of the computed first terms. If a 
numerical error had occurred in the summation for any one of the 21 ele- 
ments in each matrix, it would have definitely impaired the convergence in 
the iteration procedure. In the first eight cases of Table I, the iteration con- 
verged to within 1 % of the final value very rapidly. In the ninth case (marked 
with an asterisk in Table I), however, the iterated eigenvalue oscillated within 
+ 8% with no indication of improvements; the reason for this may have been 
some undetected numerical error which had occurred in the summation of the 
series in one or more elements of the corresponding matrix. Hence it may be 
concluded that the values of k in Table I for cases 1 to 8 are accurate within 
1% while for case 9, only within about 10%. 


d) Graphs for Computation of Buckling Loads. 

The computed values of k in Table I are plotted versus § in Fig. 2 and 
versus 2B/L in Fig. 3; the latter graph is plotted on a semi-logarithmic scale. 
Smooth curves have been drawn through the three available points for each 
curve. It can be demonstrated ! that between the three plotted points, the 
curves have to be continuous and must possess continuous derivatives; hence, 
the general shape of the curves must be that drawn in Figs. 2 and 3. 

A comparison of the curves in Fig. 2 with those applicable to uniformly 
loaded plates (Ref. (2), page 330 Fig. 172) will reveal the similarity in the 
curves for the two cases. 

One additional curve, based on the computed values of k, is drawn in 
Fig. 4. It shows the variation in strength (p,,) of a plate of constant depth 
(hy) and with constant width (2B,) of the strip loading, with the variation in 
the length (L). For the case shown in Fig. 4, the width (2B,) was arbitrarily 
so chosen that for a square plate (e.e., 8 = hy/L, = 1), 2By/Ly = 1/15. It is 
apparent that as the length (L) of the plate is varied, the ratio 2B;/L varies 
accordingly. Most of the points on the curve have been interpolated from the 
curves in Figs. 2 and 3. 
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The curve in Fig. 4 has been drawn for the instructive observation that the 
buckling strength of plates of various lengths loaded and supported as in Fig. 1, 
vary very little when they are sufficiently long, say if 8<1/3. 


e) Application to Plates With Supports Different From Figure 1. 

A common practical case is shown in Fig. 5, where the plate is supported 
at the bottom corners rather than by shearing stresses along its vertical 
edges (Fig. 1). The curves of Figs. 2 and 3 could be applied to a plate so 
supported provided the plate is long enough, say @< 1/3 (and obviously if 
buckling at the supports is prevented). This can be explained, on the one 
hand, by St. Venant’s principle, and on the other, by tests which will be dis- 
cussed in the next section. 

According to St. Venant’s principle, the manner of application of the reac- 
tions at the extreme vertical edges of the plate will not materially affect the 
stress distribution in the region of the plate sufficiently far removed from 
these edges; for practical purposes, this region can be considered to lie with- 
in the middle portion of the plate at a distance h (depth of the plate) from each 
vertical edge. The buckling load, of course, depends on the magnitude and the 
distribution of the stresses in the plate. 

In tests, failure of the plates was characterized by deformations out of the 
plane of the plates in the region directly underneath the external load. This 
is shown in Fig. 6. These deformations, in all tests, occurred in a region 
which was in the central portion of the plate; this region had a width approxi- 
mately equal to the depth of the plate, as shown in Fig. 6. 

For the above two reasons it seems justified to apply the derived curves 
to long plates supported in any manner at their vertical edges. Furthermore, 
it seems reasonable to expect that these curves would also apply to long plates 
fixed, or elastically restrained, at their vertical edges, and not necessarily 
hinged as was assumed in Fig. 1. 

The above conclusions would not hold for short plates when 8>1/3. For 
such plates, the manner in which the reactions are applied and the elastic con- 
ditions at the edges, are of paramount importance. No general solution for 
stress distribution in plates supported as in Fig. 6 has yet been derived. 
Numerical results for specific cases are available, obtained by finite differ- 
ences or similar methods.'") For these cases, however, the buckling loads 
also would have to be determined by some numerical method, for each in- 
dividual case. This was the reason for choosing a plate supported as in Fig. 1 
for which there is a known stress function (Eqn. (2)), since a general solution 
of the buckling load in terms of the plate and load parameters was desired. 

Finally it may be mentioned that the work of footnote 1, as indicated in this 
paper, can be extended to the buckling load determination of plate girder webs. 
Although some modifications in the evaluated matrices would be needed, the 
work already performed constitutes the bulk of the mathematical development 
that is required for the determination of the buckling loads for the usual cases 
of stiffened and unstiffened webs of plate girders. The theoretical procedure 
for plate girders is outlined in footnote 1. 


Ill. Tests 


a) Description of the Tests 

Tests, referred to previously in the paper, have been carried out by the 
author at Cornell University as a part of an extensive research project on 
Web Crippling of Light Gage Steel Beams for the American Iron and Steel 
Institute. A complete description of these tests has been given in reference (8), 
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A sketch of a typical test specimen and its loading is shown in Figs. 7a and 
7b. As may be seen, every test specimen has two vertical webs. If each of 
these webs is looked upon as a vertical plate, they are then in many respects 
similar to the case considered in this paper. The main differences between 
the tested plate and the theoretically treated plate are: a) in end reactions, 
and b) in shear stresses along the top and bottom edges. In the tested plate, 
due to the presence of flanges, there are horizontal shearing stresses trans- 
mitted by the flanges to the top and bottom edges which do not exist in the 
theoretically treated plate. Otherwise, the correspondence in the boundary 
conditions and in other features is very close between the two cases. This 
correspondence is discussed more fully in footnote 1. 

The test specimens were of 16 gage sheet (0.060 in. thick), 3 to 12 inches 
deep and all of constant ratio of depth to length namely, 8 = 1/4. The width 
(2B) of the central bearing plate through which the external load was applied 
(Fig. 7a) varied within the range from 2B/L = 1/65 to 1/4.8. The steel in one 
half the number of the test specimens had a yield stress of about 33000 psi 
and in the other half of about 54000 psi. 

An unavoidable phenomenon in the tests was the bending of the webs out of 
the vertical plane. This was caused by slight eccentricity of the external load 
with respect to the plane of the web due to the fact that the load was applied 
through the flanges of the specimens (Fig. 7a) which in turn had rounded cor- 
ners at junction with the web. Radius (r) of corner bents was, r = t. 

During the progress of a test at moderately high loads, but still before 
failure, the webs bent inwards (Fig. 7c) and deflected out of their plane. These 
deflections were of relatively small magnitude and extended throughout the 
depth of a web in the vicinity of the external load (i.e. central bearing plate in 
Fig. 7a). This was due, undoubtedly, to bending caused by the eccentric load- 
ing. 

At failure, there was a sudden bulging of the web with large deflections 
under, and in the immediate vicinity of, the central bearing plate, as shown in 
Fig. 6. 

It will be realized from the above that although the webs were subjected to 
bending, buckling was the governing factor in failure. Nevertheless, the 
presence of bending is bound to have its effect and should be borne in mind in 
comparing test failure loads with the theoretical buckling loads. Because of 
the effect of bending, one would expect higher failure loads for the higher 
strength steel than for the lower strength steel. Due to bending, the absolute 
depth (h) of the web would also have an additional effect on the failure load. 
Finally, the membrane action of the web (post buckling strength) should not be 
neglected in comparing failure loads with the theoretical buckling loads. 


b) Test Results; Comparison With Theoretical Buckling Loads. 

The failure loads obtained in the tests are plotted in Fig. 8 to the same 
scale as the theoretical curve for buckling loads. The latter is shown by a 
solid line and was obtained from Fig. 3. The broken lines join tested points 
of identical test specimens but of varying 2B/L. 

Inspection of Fig. 8 will reveal the following: 


1) The trend of the test failure loads is the same as that of the theoretical 
curve. 

2) Agreement between the tests and the theoretical loads is closer for the 

specimen of the lower yield steel than for the higher yield steel. The 

agreement is closest for the 6" and the 9" deep specimens (h/t of 90 

and 150, respectively) of the lower yield group specimens. 
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3) The theoretical curve predicts mostly buckling loads smaller than test 
failure loads. This is particularly true for deeper specimens of high 
yielding steel. 


The discrepancies between the theoretical curve and the tests could be at- 
tributed to the following: 


a) Difference in the experimental set-up and the theoretical case analysed. 

b) Post-buckling strength. 

c) Use of only six terms in Eqn. (10). Because buckling in the treated case 
is mostly a localised effect, greater accuracy should be expected by us- 
ing more than six terms in the expansion of the expression for the 
buckled shape. However, addition of one or more terms to Eqn. (10) 

would increase considerably the mathematical work involved. 


CONCLUSION 


Formula for the buckling load of flat rectangular plates loaded with strip 
loading at one edge has been derived. This formula has been reduced to the 
familiar form for the buckling load of flat rectangular plates loaded uniformly 
along the two opposite edges; the difference between the formulae for the two 
cases is in the buckling coefficient k. In the formula derived in this paper, 

k depends on the plate dimensions and loading widths. Curves for the deter- 
mination of this coefficient have been plotted and can be readily used for the 
design purposes. 
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8. L. Zetlin and G. Winter: 65th and 66th Progress Reports to the American 
Iron and Steel Institute on Tests on Light Gage Steel Beams of Cold Formed 
Steel, Cornell Univ., Sept. 1952 and Jan. 1953. 
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PROCEEDINGS PAPERS 


The technical papers published in the past year are presented below. Technical-division 
sponsorship is indicated by an abbreviation at the end of each Paper Number, the symbols 
referring to: Air Trarsport (AT), City Planning (CP), Construction (CO), Engineering Mechanics 
(EM), Highway (HW), Hydraulics (HY), Irrigation and Drainage (IR), Power (PO), Sanitary 
Engineering (SA), Soil Mechanics and Foundations (SM), Structural (ST), Surveying and Mapping 
(SU), and Waterways (WW) divisions. For titles and order coupons, refer to the appropriate 
issue of “Civil Engineering” or write for a cumulative price list. 


VOLUME 80 (1954) 


SEPTEMBER: 484(ST), 485(ST), 486(ST), 487(CP)°, 488(ST)°, 489(HY), 490(HY), 491(HY)°, 
492(SA), 493(SA), 494(SA), 495(SA), 496(SA), 497(SA), 498(SA), 499(HW), 500(HW), 501(HW)°, 
502(WW), 503(WW), 504(Ww)°, 505(CO), 506(CO)©, 507(CP), 508(CP), 509(CP), 510(CP), 
511(CP). 


OCTOBER: 512(SM), 513(SM), 514(SM), 515(SM), 516(SM), 517(PO), 518(SM)°, 519(IR), 520(R), 
521(IR), 522(IR)°, 523(AT)°, 524(SU), 525(SU)°, 526(EM), 527(EM), 528(EM), 529(EM), 
530(EM)°, 531(EM), 532(EM)-, 533(PO). 


NOVEMBER: 534(HY), 535(HY), 536(HY), 537(HY), 538(HY)°, 539(ST), 540(ST), 541(ST), 542(ST), 
543(ST), 544(ST), 545(SA), 546(SA), 547(SA), 548(SM), 549(SM), 550(SM), 551(SM), 552(SA), 
553(SM)°, 554(SA), 555(SA), 556(SA), 557(SA). 


DECEMBER: 558(ST), 559(ST), 560(ST), 561(ST), 562(ST), 563(ST)°, 564(HY), 565(HY), 566(HY), 
567(HY), 568(HY)°, 569(SM), 570(SM), 571(SM), 572(SM)°, 573(SM)°, 574(SU), 575(SU), 576(SU), 
577(SU), 578(HY), 579(ST), 580(SU), 581(SU), 582(Index). 


VOLUME 81 (1955) 


JANUARY: 583(ST), 584(ST), 585(ST), 586(ST), 587(ST), 588(ST), 589(ST)°, 590(SA), 591(SA), 
592(SA), 593(SA), 594(SA), 595(SA)~, 596(HW), 597(HW), 598(HW)°,599(CP), 600(CP), 601(CP), 
602(CP), 603(CP), 604(EM), 605(EM), 606(EM)°, 607(EM). 


FEBRUARY: 608(WW), 609(WW), 610(WW), 611(WW), 612(WW), 613(WW), 614(WW), 615(WW), 
616(WW), 617(1R), 618(1R), 619(1R), 620(IR), 621(1R)°, 622(1R), 623(1R), 624(HY)©, 625(HY), 
626(HY), 627(HY), 628(HY), 629(HY), 630(HY), 631(HY), 632(CO), 633(CO). 


MARCH: 634(PO), 635(PO), 636(PO), 637(PO), 638(PO), 639(PO), 640(PO), 641(PO)*, 642(SA), 
643(SA), 644(SA), 645(SA), 646(SA), 647(SA)°, 648(ST), 649(ST), 650(ST), 651(ST), 652(ST), 
653(ST), 654(ST)°, 655(SA), 656(SM)°, 657(SM)°, 658(SM)°. 


APRIL: 659(ST), 660(ST), 661(ST)°, 662(ST), 663(ST), 664(ST)°, 665(HY)°, 666(HY), 667(HY), 
ener 669(HY), 670(EM), 671(EM), 672(EM), 673(EM), 674(EM), 675(EM), 676(EM), 677(EM), 


MAY: 679(ST), 680(ST), 681(ST), 682(ST)©, 683(ST), 684(ST), 685(SA), 686(SA), 687(SA), 688(SA), 
689(SA)©, 690(EM), 691(EM), 692(EM), 693(EM), 694(EM), 695(EM), 696(PO), 697(PO), 698(SA), 
699(PO)°, 700(PO), 701(ST)°. 


JUNE: 702(HW), 703(HW), 704(HW)°, 705(IR), 706(IR), 707(IR), 708(IR), 709(HY)°, 710(CP), 
711(CP), 712(CP), 713(CP)®, 714(HY), 715(HY), 716(HY), 717(HY), 718(SM)°, 719(HY)¢, 
720(AT), 721(AT), 722(SU), 723(WW), 724(Ww), 725(WW), 726(Ww)°, 727(WW), 728(IR), 
729(IR), 730(SU)°, 731(SU). 


JULY: 1732(ST), 733(3T), 734(ST), 735(ST), 736(ST), 737(PO), 738(PO), 739(PO), 740(PO), 
741(PO), 742(PO), 743(HY), 744(HY), 745(HY), 746(HY), 747(HY), 748(HY)°, 749(SA), 750(SA), 
TOL(RA), 753(SM), 754(SM), 755(SM), 756(SM), 757(SM), 758(CO)°, 759(SM)°, 


AUGUST: 761(BD), 762(ST), 763(ST), 764(ST), 765(ST)©, 766(CP), 767(CP), 768(CP), 769(CP), 
TT0(CP), 771(EM), 772(EM), 773(SA), 774(EM), 775(EM), 776(EM)°, 777(AT), 778(AT), 
779(SA), 780(SA), 781(SA), 782(SA)°, 783(HW), 784(HW), 785(CP), 786(ST). 


SEPTEMBER: 787(PO), 788(IR), 789(HY), 790(HY), 791(HY), 792(HY), 793(HY), 794(HY)‘, 
795(EM), 796(EM), 797(EM), 798(EM), 799(EM)°, 800(WW), 801(WwW), 802(WW), 803(WwW), 
804(WW), 805(WW), 806(HY), 807(PO)°, 808(IR)°. 


c. Discussion of several papers, grouped by Divisions. 
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